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Radio continuum and CO emission in star-forming galaxies. 
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Abstract. We combine the radio continuum images from the NRAO VLA Sky Survey with the CO-line observa- 
tions from the extragalactic CO survey of the Five College Radio Astronomy Observatory to study the relationship 
between molecular gas and the star formation rate within the disks of 180 spiral galaxies at 45" resolution. We 
find a tight correlation between these quantities. On average, the ratio between the radio continuum and the CO 
emission is constant, within a factor of 3, both inside the same galaxy and from galaxy to galaxy. The mean star 
formation efficiency deduced from the radio continuum corresponds to convert 3.5% of the available molecular 
gas into stars on a time scale of 10 8 yr and depends weakly on general galaxy properties, such as Hubble type 
or nuclear activity. A comparison is made with another similar analysis performed using the H Q luminosity as 
star formation indicator. The overall agreement we find between the two studies reinforces the use of the radio 
luminosity as star formation rate indicator not only on global but also on local scales. 
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1. Introduction 

One of the major goals of the studies of external galaxies 
is understanding the relationship between the star forma- 
tion rate (SFR) and the physical condition in the inter- 
stellar medium. Several indicators have been suggested to 
estimate the SFR (of massive stars) in galaxies. These 
include the U-band magnitude, the strength of Balmer 
lines emission, the far-infrared (FIR) emission and the ra- 
dio luminosity; the rates inferred from the different indi- 
cators span almost four orders of magnitude going from 
1CT 2 to 10 2 M yr" 1 . Cram et al. (1998) checked the 
consistency between the SFR deduced from the U-band, 
H Q , FIR and radio luminosity using a sample of 700 local 
galaxies. They noted that there are systematic differences 
between these various indicators. In particular they sug- 
gested that the H Q luminosity may underestimate the SFR 
by approximately an order of magnitude when the SFR 
is > 20 Mq yr -1 . They concluded that the radio contin- 
uum luminosity at decimeter wavelengths of a star forming 
galaxy provides a better way to estimate the current rate 
of star formation. The radio continuum emission at 1.4 
GHz from a star-forming galaxy is mainly synchrotron ra- 
diation produced from relativistic electrons accelerated by 
supernovae explosions (Lequeux 1971). Indeed the radio 



Send offprint requests 
murgia@ira . bo . cnr . it 



to: 



M. 



Mursia, 



continuum luminosity appears to be directly proportional 
to the rate of formation of supernovae (Condon 1992). 
This view is reinforced by the tight correlation existing 
between the radio luminosity and the FIR for spiral galax- 
ies (see e.g. Condon 1992). Since the radio continuum at 
1.4 GHz does not suffer significant extinction, the radio 
luminosity constitutes a very useful tool to determine the 
current SFR in a spiral galaxy. 

Since the discovery that stars form in molecular clouds, 
it is essential to determine, not only the rate, but also the 
efficiency of conversion of the interstellar gas in stars; i.e. 
the star formation efficiency (SFE). The SFE measures 
the formation rate of young stars per unit of mass of gas 
available to form those stars. Determining the SFE is im- 
portant to distinguish a situation in which a high SFR in- 
dicates a higher efficiency in converting gas in stars rather 
than a higher gas quantity. 

The CO molecule luminosity and the virial mass of 
giant molecular clouds correlate very well in our Galaxy 
and in other nearby spirals (Young & Scoville 1991 and 
references therein) . The comparison of different SFR trac- 
ers with the mass of molecular clouds provides indeed an 
important tool to investigate the behaviour of the SFE 
within and among galaxies. Many studies have been con- 
cerned with the behaviour of the star formation process on 
global scales, averaged over the entire star-forming disk. 
These works showed that the disk-averaged star forma- 



2 



M. Murgia et al.: Radio continuum and CO emission in star-forming galaxies. 



tion process is well described by a Schmidt (1959) law 
of the type Ssfr oc S^ s , where Esfr and S gas are the 
observable surface density of SFR and total (atomic + 
molecular) gas density, respectively, and the exponent N 
typically ranges from 1.3 to 1.5 (Kennicutt 1998). 

An interesting development of these global studies, 
the investigation of the behaviour of the SFE within the 
disks of the individual galaxies, provides much physical 
insight into the star formation process itself. The extra- 
galactic CO survey of the Five College Radio Astronomy 
Observatory (Young ct al. 1995, hereafter FCRAO CO 
Survey) provided a uniform database of CO data for 300 
galaxies at a resolution of 45", opening the possibility to 
extend the study of the Schmidt relationship of the SFR 
versus the H 2 density over the same physical regions well 
inside the galaxy disks. Since the star formation process 
involves the molecular gas directly, some authors recog- 
nized that the determination of the Schmidt law assumes 
a clear physical meaning if restricted to this gas compo- 
nent. Moreover, in the considered regions the molecular 
gas is dominant over the atomic one and, contrary to this 
latter, its azimuthally averaged distribution follows closely 
the radial profiles of the main SFR indicators (Tacconi 
and Young 1986, Young & Scoville 1991). Rownd & Young 
(1999; hereafter RY99) conducted an H a imaging of 121 of 
these galaxies, determining the local relationship between 
the SFR and the molecular gas. They found a correla- 
tion between these two quantities and concluded that for 
face-on spiral, in general, there are no strong SFE gradi- 
ents across the star-forming disks. The majority of large 
SFE variations they found are seen between adjacent disk 
points, reflecting regional differences in the SFE, and any 
radial gradients are at most a secondary effect. In con- 
trast, they pointed out that consistent radial variations 
(up to an order of magnitude or more) of the SFE exist 
within many highly inclined galaxy disks. They attributed 
the decreasing SFE towards the centers of these galaxies 
to a large amount of dust extinction on the H Q luminosity. 

Adler ct al. (1991) found a correlation between the 
radio continuum flux density at 20 cm and the CO line 
emission on global scales for a sample of 31 spiral galaxies. 
They also studied the relationship of these two quantities 
within the disks of 8 nearby well resolved spiral galaxies, 
finding that their ratio is constant both inside the same 
galaxy and from galaxy to galaxy. 

The work we present here is complementary to the 
analysis of RY99 and extends that of Adler et al. (1991). 
We combined the radio continuum images at 1.4 GHz 
from the NRAO VLA Sky Survey (NVSS, Condon ct al. 
1998) with the FCRAO CO survey to study the rela- 
tionship between the radio continuum and the molecular 
gas point-to-point within the disks of 180 star-forming 
spiral galaxies. It is important to stress that we are 
comparing two homogeneous data set with the same 
angular resolution of 45". 



The paper is organized as follows: in Sect. 2 and 
Sect. 3 we present the sample used and we describe the 
data analysis, respectively. In Sect. 4 we present the 
results of the statistical analysis and in Sect. 5 we discuss 
the results obtained. 

We use a Hubble constant Ho=50 km s _1 Mpc _1 
throughout the paper. 

2. Sample selection 

To investigate the star formation law within the disks of 
normal galaxies, we combined the data from two public 
surveys. We use the NVSS for the radio continuum and 
the FCRAO CO survey for the molecular gas, respectively. 

The NVSS was performed at 1.4 GHz with the Very 
Large Array (VLA) in D configuration. It has an angular 
resolution of 45" (FWHM), a noise level of 0.45 mJy/beam 
(la) and covers all the sky north of declination -40°. The 
shortest baseline is 35 m, corresponding to ~167A, there- 
fore structures up to about 10' in angular size are properly 
imaged. 

The FCRAO CO survey comprises 300 galaxies ob- 
served along the major axis of the disk for a total of 1412 
locations. Most of the galaxies in the survey are spirals or 
irregulars north of declination -25°. At the frequency of 
the CO J = 1 - transition (115.27 GHz) the FWHM 
of the 14-m FCRAO telescope is 45". The weakest line 
detected depends on the width of the line, and hence on 
the velocity field within the beam. The uncertainties on 
the individual line intensity vary from galaxy to galaxy. A 
conservative estimate of the rms noise, including the cal- 
ibration, baseline removal, and the rms noise per channel 
is about 25% (the median signal-to-noise ratio is 4). 

We note that the two surveys have uniform sensitiv- 
ity and identical angular resolution. This fact circumvents 
the difficulties deriving from the comparison of data from 
multiple instruments or studies which are subtly incom- 
patible either because of inconsistent signal-to-noise ratios 
or unmatched resolution. 

The original FCRAO CO survey includes 300 galax- 
ies selected from the RC2 (de Vaucouleurs et al. 1976) or 
the IRAS database satisfying at least one of the following 
criteria: i) < 13.0, ii) Sao > 5 Jy or iii) S*ioo > 10 
Jy. Although the FCRAO CO survey is not a complete 
sample in terms of flux-density or volume limit, the ob- 
served galaxies cover a wide range of luminosity, morphol- 
ogy and environments. For this reason they represent an 
ideal database to study the behaviour of the star forma- 
tion process and the molecular gas in a wide variety of 
conditions. 

Since our interest was primarily to investigate the be- 
haviour of star formation within the galaxy disks, we have 
selected, from the FCRAO CO survey, a sub-sample of 180 
objects for which there were at least three different obser- 
vations of the CO line in the disk. 

Distributions of morphological types, angular diame- 
ters, distances and linear diameters of the sample along 
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Fig. 1. Distributions of morphological types (top left panel), optical angular diameters (top right panel), distances 
(bottom left panel) and linear diameters (bottom right panel): the solid portion of the histograms indicates the sub- 
sample used in this work (180 galaxies) whilst the empty portion indicates the distribution of the FCRAO CO survey 
(300 galaxies). Because of its outstanding angular size, the galaxy NGC598 (M33, D25 = 62') is not included in the 
histogram of angular diameters. 



with the corresponding distributions for the FCRAO CO 
survey are shown in Fig. ^. 

The galaxies in our sample have morphological types 
ranging between SO and 10. Most of them (89%) are spiral 
galaxies with morphological type between Sa-Sd. 

The majority of the galaxies (about 90%) have an op- 
tical angular diameter (from RC2) D25 <10', the median 



angular diameter being about 5'. This ensures that for 
all of them the NVSS properly recovered the flux density 
of the extended structures (see above) . The remaining 20 
galaxies have an optical angular diameter from 10' to 25'. 
For these galaxies it is possible that the NVSS missed a 
significant fraction of flux from the extended structure. 
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Because of its outstanding angular size, NGC598 (M33, 
Z?25 = 62') had been excluded from the analysis. 

The distances of the galaxies in our sample (taken from 
Young et al. 1995) span from the Local Group up to about 
80 Mpc. Over 87 galaxies are at the distance of the Virgo 
Cluster (20 Mpc). 

The linear diameters range from ~ 4 to ~ 100 kpc. 
The median value is 31 kpc. 

Our selection excluded most of the galaxies with an 
angular diameter less than 3' , i.e. the intrinsically small 
galaxies (linear diameter smaller than 30 kpc) and the 
more distant ones. 

The complete list of the galaxies in our study 
(including the galaxy name, Hubble type, incli- 
nation, angular diameter and distance) is avail- 
able in electronic format at http://www.ira.bo.cnr.it/ 



~crapsi_s/RADIOCO/Article/tab_art.txt. 
3. Data analysis 

For each galaxy in our sample we extracted the radio sur- 
face brightness, B1.4, from the NVSS images . We averaged 
the NVSS pixel values within circular regions of 45" di- 
ameter, centered on the same pointing positions of the 
FCRAO survey. Following the same formalism adopted 
by RY99 for the H a emission, the radio brightness is ex- 
pressed in terms of the ratio between the radio luminosity 
at 1.4 GHz, L1.4, and the de-projected surface, E, inter- 
cepted by the beam on the galaxy disk: 



L1.4 (W Hz" 1 ) 



3.2 x 10 18 B 1A 



( mJy \ 

\beam J 



cos 1 



(1) 



E (kpc 2 ) 

where i is the inclination angle of the galaxy with respect 
to the plane of the sky. The surface luminosity has the 
advantage to be distance independent as opposed to the 
global luminosity. 

Since most of continuum radio luminosity of normal 
galaxies is produced by relativistic electrons accelerated 
by supernovae explosion and the supernova rate (z^sn) is 
directly related to the SFR of massive stars, a relation is 
expected between the radio luminosity and the star for- 
mation rate. In the following we indicate with SFR the 
formation rate of stars with mass M > 5M . Condon 
(1992) calibrated empirically the ^sn-SFR relation using 
the supernova rate and the radio luminosity of our Galaxy: 

L v (W Hz" 1 ) = 5.3 x 10 2 V G g z • SFR (Moyr^ 1 ) (2) 

where a Miller-Scalo (1979) initial mass function (IMF) 
is assumed and a is the radio spectral index (S v cx v~ a ). 
Condon (1992) noted that this relation should apply to the 
majority of normal galaxies in order to be consistent with 
the tight FIR/radio correlation widespreadly observed. 

Assuming for the radio spectral index a typical value 
of a — 0.8, from Eq. (|l|) and Eq. (||) the relation between 
the star formation rate per unit surface, Esfr, and the 
radio brightness at 1.4 GHz is found to be: 



Mr. 



-'SFR 



v yr kpc 2 



3.0 x 10- 4 Bl 4 ( £*L) 
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cos 1 
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In order to derive the H2 surface density from the 
FCRAO CO survey integrated CO intensity Iqo , we used 
the formula (RY99): 



E H2 (M pc" 2 ) = 8.5 Jco(K km s" 1 ) cos 



(4) 



In this formula it is assumed a constant linear conversion 
factor between the CO luminosity and H2 mass equal to 
V co = 2.8 x 10 20 H 2 cnr 2 [K km s" 1 ]" 1 (Bloemen et al. 
1986). 

A realistic estimate of the uncertainties in both Esfr 
and Eh 2 surface densities should consider several system- 
atic effects. The relation between the SFR and the radio 
luminosity is based on many not well proved assumptions, 
such as the IMF thresholds and slope and the extrapola- 
tion of the ^sn-SFR Milky Way relation to other galaxies. 
Cram et al. (1999) pointed out that different modelling of 
these parameters introduce scaling uncertainties up to a 
factor of 2. The dominant errors in the gas density are the 
variation on the CO/H2 conversion factor. These varia- 
tions can be as high as ±40% for luminous spiral galaxies 
as those studied in this work (Devereux & Young 1991). 
Despite these uncertainties, the data provide very strong 
constrains on the form of the SFE because of the wide 
ranges of SFR and gas densities explored. 

Consistently with the definition given by RY99, the 
radio SFE is calculated by the ratio of Esfr to Eh 2 ■ 
In terms of our observables (-B1.4 and Iqo ) the SFE is 
expressed by 



SFE (yr- 1 ) 



^SFR 

Eh 2 



9.4 x 10 



-11 B 1A 
Ico 



(5) 



where Si. 4 and Ico are measured in mJy/beam and 
K km s _1 , respectively. 

The SFE defined by Eq. (g) gives the fraction of molec- 
ular gas converted to massive stars per year. Since the 
typical lifetime of the synchrotron radiating electrons is 
shorter than 10 8 yr (Condon 1992), the SFR inferred 
from the radio luminosity traces a stellar population not 
older than this timescale (hereafter t$). The percentage of 
molecular gas consumed over all this period is 



SFE rs (%) = 0.94 



4. Results 



Ico 



(6) 



The data analysis presented in Sec. 3 allows one to com- 
pare the SFR and the molecular gas content to determine 
the star formation efficiency along the galaxy disks. 

Overlays of optical (grayscale) and radio continuum 
(contours) images for eight selected galaxies^] are shown 
in the left columns of Figs. || and || where the circles 



The plots for all the galaxies in the sample can be found 

at the web page: 

|http: //www. ira.bo. cnr.it/~crapsi_s/RADIOCO/Catalog 
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Fig. 2. Examples of selected galaxies, see text. Solid and short-dashed lines in middle panels refer to the mean and 
standard deviation of the SFE Tg of the entire sample. Solid and dashed lines in right panels refer to Ico an Bia 
profiles, respectively. 
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Fig. 3. Examples of selected galaxies, see text. Solid and short-dashed lines in middle panels refer to the mean and 
standard deviation of the SFE Tg of the entire sample. Solid and dashed lines in right panels refer to Ico an Bia 
profiles, respectively. 
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Fig. 4. Histogram of star formation efficiency at 628 
pointing positions for the entire sample. Only points with 
a signal-to-noise ratio greater than 2a have been consid- 
ered. Mean and dispersion refer to the logarithmic val- 
ues of SFE rs , see text. Values greater than 100% indicate 
present SFRs that cannot be supported for more than 10 8 

yr- 



indicate the positions and the beam size of the CO obser- 
vation. In the middle columns we plot B\ a versus Ico re- 
porting also the corresponding values of Ssfr and £h 2 in 
the upper and right axis, respectively. The panels in the 
right columns show the radio brightness (Si. 4), the CO 
integrated intensity (Ico ) an d their ratio as a function 
of distance from the galaxy center. The convention is that 
radius is positive for positive right ascension (or declina- 
tion) pointing shifts. The optical images are taken from 
the red Palomar Digitized Sky Survey. The NVSS radio 
contours start at 0.9 mJy beam^ 1 (2a) and are spaced 
by a factor of y/2. In the plots, error bars and arrows in- 
dicate respectively measurement uncertainties and upper 
limits (2a). In the middle column panels, the reference 
lines represent the mean (solid) and standard deviation 
(short-dashed) of the SFE computed using all the detec- 
tions in the whole sample (see Sect. 4.2). In right column 
panels, the dashed and continuous lines show the -B1.4 and 
Ico trends, respectively. 

4.1. The SFE within galaxy disks 

Most of the galaxies in our sample are at distances between 
10 and 40 Mpc, the corresponding linear resolution of the 
45" beam is between 2.2 and 8.7 kpc, i.e. so large to sample 
the contribution from many molecular clouds complexes. 
In the following, one should keep in mind that the star- 



forming regions we are observing are comparable with the 
width of the spiral arms only for the nearest objects. 

The galaxies shown in Figs. 2 and 3 are representative 
of the diversity of behaviour seen in the distributions of 
-B1.4 and Ico within galaxy disks. 

The most striking feature is the linear correlation be- 
tween these two quantities observed for many edge-on 
and face-on galaxies (see Fig. 2). In these cases, -B1.4 and 
Ico present the same scaling from the galaxy center out- 
ward, resulting in a constancy of the SFE along the disk. 
However, there are clear examples of disks characterized 
by systematic SFE trends (see Fig. 3). A _Bi.4 gradient 
steeper (flatter) with respect to Ico one implies a SFE 
decreasing (increasing) with radius (e.g. NGC5236 and 
NGC5247). Calculating the SFR from non-thermal radio 
continuum allows us to include in the analysis high in- 
clined galaxies, such as NGC1055 and NGC3079, which 
generally suffer from extinction in the optical band (see 
Sect. 5.1). By fitting a power law of the form B\ A cx Ico N , 
we found that the fraction of linear (0.5 < N < 2) corre- 
lations is 67%, while the fractions of super-linear (N > 2) 
and sub-linear (N < 0.5) correlations are 23% and 10%, 
respectively. We examine the composite correlation in- 
cluding all the pointings in the sample in Sect. 4.3. 

4.2. The SFE variation among galaxies 

Having established a correlation between the radio contin- 
uum surface brightness and the CO line intensity within 
galaxies, we investigated the variation of the SFE among 
galaxies. Fig. 4 shows the distribution of the SFE T8 for the 
galaxy sample. SFE Ts varies from about 0.1 up to more 
than 100%. Values greater than 100% indicate present 
SFRs that cannot be supported for more than 10 8 yr. Since 
the distribution of logarithmic values is more Gaussian we 
calculated the mean and the dispersion of log(SFE rs ). The 
mean of log(SFE T8 ) is approximately 3.5% of the gas con- 
sumed per 10 8 yr and the rms of the distribution is slightly 
less than a factor of 3. SFE Ts mean and dispersion of the 
sample are traced in middle plots of Figs. 2 and 3 as solid 
and dashed reference lines, respectively. 

We examined also the variation of SFE T8 among galax- 
ies compared to the morphological type (see Fig. 5). The 
mean star formation efficiency varies weakly (about 25%) 
with the morphological type going from SO to Scd galaxies. 

Larger SFE close to the galaxy centers might be at- 
tributed to the presence of an active nucleus (AGN). 
Fig. 6 shows the distribution of disk and nuclear SFE T8 in 
Seyfert and normal galaxies. Seyferts have a slightly higher 
SFE Ts than normal galaxies. The mean of the log SFE T8 , 
for Seyfert and normal galaxies, is: 4.1% and 3.3% in 
the disks and 5.5% and 4.3% in the nuclei, respectively. 
Most Seyferts have a nuclear SFE T8 comparable with 
normal galaxies. Only in few cases the nuclear emis- 
sion is dominated by the radio source related with the 
AGN, e.g. NGC1068 (Wilson & Ulvestad 1987), NGC4151 
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Fig. 5. Histograms of star formation efficiency sorted 
by galaxy type. Only points with a signal-to-noise ratio 
greater than 2a have been considered. Mean and disper- 
sion refer to the logarithmic values of SFE Tg , see text. 
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Fig. 6. Histograms of star formation efficiency for nu- 
clear (top panel) and disk (bottom panel) pointings. The 
solid line and solid portion of the histograms indicate 
normal and Seyfert galaxies, respectively. Three galaxies 
have nuclear SFE Tg greater than 100%; namely NGC2655, 
NGC3310 and NGC4151. NGC2655 and NGC4151 are 
classified as Sy2 and Syl.5 and have nuclear SFE Ts of 
100% and 170%, respectively. The galaxy NGC3310 is a 
peculiar starburst galaxy characterized by an exceptional 
nuclear SFE Tg of about 135% 



(Pedlar et al. 1993), NGC2655 (Keel & Hummel 1988) and 
NGC4388 (Irwin et al. 2000). We conclude that the AGN- 
related emission affects only marginally the estimate of 
the SFE Tg in the nuclear pointings for most Seyfert in our 
sample. 

We further investigate the behaviour of the SFE with 
respect to the galaxy inclination and size, and beam lin- 
ear resolution. Fig. 7 shows the maximum variation of 
the SFE Tg , defined as SFE max /SFE min , inside each galaxy. 
Most galaxies present SFE variation up to a factor 6, the 
median variation being a factor of 2.5. 

Fig. 7 shows that the internal SFE variations are not 
strong function of the galaxy inclination or linear resolu- 
tion of the observing beam. Considerable SFE variations 
occur also for size of the beam greater than 5 kpc, i.e. 
smoothing over large regions of the galaxy disks. This fact 
was already noted by RY99. Ten galaxies show an inter- 
nal SFE variation greater than a factor of 10. These are: 
the circumnuclear starbursts IC342, NGC253, NGC520, 
NGC660, NGC2146 and NGC3034 (see Kennicutt 1998); 
the Seyfert galaxies NGC2841 and NGC3368; the pecu- 
liar galaxy NGC3628; the HII galaxy NGC6503. M82, 
NGC520 and NGC6503, show exceptional internal SFE 



variations larger than a factor of 30. In particular, the 
nearby starburst M82 show a variation of about 2.6 order 
of magnitude. 

Finally, we investigated the behaviour of the star for- 
mation efficiency as a function of the distance from the 
galaxy centers. Fig. 8 shows the SFE as a function of ra- 
dius for all galaxies. The SFE is found to be approximately 
constant at all radii. 

4.3. The composite radio Schmidt law 

Our data offer the possibility to determine the spatially 
resolved Schmidt law using the radio continuum as SFR 
indicator for a sample of galaxies six times larger than the 
one of Adler et al. (1991). Fig. 9 shows the relationship 
between the SFR and the H2 surface densities for all the 
pointings of the sample. Dashed and dotted reference lines 
indicate the mean and standard deviation of the SFE rg for 
all the detection in the sample (see Sect. 4.1). The galax- 
ies NGC4302, NGC4526, NGC4710 and NGC4866 have 
been omitted since they have inclination i — 90°. A clear 
correlation is evident over more than three orders of mag- 
nitude, albeit a consistent scatter of 0.5 dex is present. We 
stress once more that, since this is a correlation between 
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two brightness, we can rule out any subtle effect intro- 
duced by the distance. This result confirms the correlation 
found by Adler et al. (1991) on global scales. The best 
weighted fit of the radio Schmidt law, Esfr — ax Eh^, 
is: a = 2.6(±1.2) x 1(T 4 and N = 1.3 ± 0.1. The value 
for the exponent of the radio Schmidt law is fully consis- 
tent with both the global (e.g. Kennicutt 1998) and local 
(RY99) slopes found using the H a line emission as SFR 
indicator. This exponent for the Schmidt law indicates a 
SFE which weakly increases with the SFR. 



5. Discussion and conclusions 

The most striking result of our study is the overall consis- 
tency of the SFE obtained from the 20 cm radio contin- 
uum as SFR indicator instead of the H a emission (RY99). 
Basically we found all the essential features observed by 
RY99: the constancy of the SFE both within and among 
galaxies; the weak dependence of the SFE on morpho- 
logical type and linear size of the observing beam; the ex- 
tent and slope of the composite Schmidt law. Nevertheless, 
some important differences arise. In this section we present 
in details the differences between the two works, discussing 
their nature and implications. 
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Fig. 8. Star formation efficiency as a function of radius 
for all galaxies. Limits are not shown. All the detections 
are represented as small points. Bold points and error bars 
represent the median value and the 50% of objects in each 
bin, respectively. 



5.1. Star formation rate surface densities 
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Fig. 7. Star formation efficiency variations within each 
galaxy as a function of inclination (top panel) and lin- 
ear resolution of the 45" beam (bottom panel). Seyfert 
galaxies are distinguished with open dots. The majority 
of galaxies have internal SFE variation less than a factor 
of 3. The starburst galaxies NGC3034 (M82) and NGC520 
and the HII galaxy NGC6503, have exceptional internal 
SFE variations (larger than a factor of 30). 



Cram et al. (1998) checked the consistency between the 
SFR deduced from the radio luminosity with the rates pre- 
dicted by other indicators, in particular by the H Q lumi- 
nosity. They concluded that the rates deduced from the ra- 
dio continuum and the H a luminosities, although in broad 
agreement, are affected by two systematic deviations: re- 
spectively at low (< 0.1 Moyr -1 ) and high (> 20 Moyr" 1 ) 
star formation rates, the H Q luminosity overestimate and 
underestimate the SFR deduced by the radio luminosity. 
The deviation at low rates can be attributed in part to 
problems related to the zero-point H Q luminosity calibra- 
tion. They suggest that the deviation at high SFR could 
be attributed to a larger amount of extinction by dust in 
those objects undergoing strong star formation or to par- 
ticular IMFs that favour low mass supernovae progenitors 
rather than high mass stars responsible for the H Q . 

By comparing our data set with that of RY99 we have 
the possibility to extend the consistency check between 
the surface star formation rate densities deduced from the 
radio continuum and H Q emission over the same regions 
of galaxy disks. 

Using Eq. (2) of Cram et al. (1998), we calculated the 
SFR surface density from the H Q brightness reported by 
RY99 through the formula: 



EsFR 



Mr; 



yr kpc 2 



2.6 x 10~ 2 /i(H Q ) 



(7) 



where the (de-projected) H Q surface brightness, /i(H a ), is 
measured in L Q pc~ 2 . Eq. (2) of Cram et al. (1998) ac- 
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Fig. 9. Spatially resolved radio Schmidt law for all the 
detections in the sample. Arrows indicate upper limits at 
2cr. The solid line is the best fit of the composite radio 
Schmidt law: £sfr oc £h 2 ' ■ Dashed and dotted lines re- 
fer to the mean and standard deviation of the SFE rg for all 
detections in the sample: 3.5% and a factor of 2.8 respec- 
tively. The galaxies NGC4302, NGC4526, NGC4710 and 
NGC4866 have been omitted since they have inclination 
i = 90°. 



counts for a correction of 1.1 magnitude in the luminosity 
to compensate for the mean extinction in H Q (Kennicutt 
1983a). Fig. 10 shows the direct comparison between the 
radio and H Q SFRs inferred respectively from Eq. (U) and 
Eq. (0) for 453 position within 102 galaxies of the sam- 
ple. All galaxies are shown in the top panel of Fig. 10. 
In general there is a broad agreement between the two 
SFR indicators. Considering the distribution of the log- 
arithm values of the ratio R = Esfr(h q )/S S fr(i.4GHz), 
we have that R have a mean of 0.83 with a standard de- 
viation of 2.8, i.e. the H a is underestimating the surface 
SFR as compared to continuum radio emission. Sorting 
the sample by galaxy inclination indicates that this ef- 
fect is strongly correlated with galaxy orientation. For the 
face-on subsample (i < 40°) R has a mean of 1.1 and 
a dispersion of a factor 2.3. At intermediate inclinations 
(40° < i < 70°) mean and dispersion of R are 0.8 and a 
factor of 3, respectively. Finally, for highly inclined galax- 
ies (i > 70°) mean and dispersion of R are 0.7 and a factor 
of 2.9, respectively. In particular, SFR surface densities 
deduced by the H Q for highly inclined galaxies are sys- 
tematically underestimated for Xsfr(i.4 GHz) greater than 
~ 3 x 10~ 3 M Q yr -1 kpc -2 . Fig. 11 shows two edge-on 
spirals, NGC3079 and NGC1055 (see Figs. 2 and 3) and 
two well known starburst galaxies, NGC2146 and M82, 
in which SFRs deduced by the H a emission are underes- 
timated by more than a factor of about 10. As a conse- 
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Fig. 10. Comparison of SFRs surface densities deduced 
from 1.4 GHz luminosity (horizontal axis) and H Q emis- 
sion (vertical axis) for 102 galaxies of our sample in com- 
mon with RY99. The various panels show different val- 
ues of galaxy inclinations i, as indicated. There are 453 
detections, 65 of which are upper limits (at 2a). The 
reference dashed lines indicate 1:1 relations. The sub- 
sample of face-on galaxies (i < 40°) shows the best 
correlation, with a dispersion of a factor of 2.3. SFRs 
surface densities deduced by the H Q for highly inclined 
galaxies (i > 70°) are systematically underestimated for 
EsfrI-4 GHz > 3 x IO" 3 M Q yr^kpc" 2 . 



quence, the SFE variations deduced for these objects by 
RY99, on the basis of H Q observations, can be attributed 
to extinction, as suspected by these authors. 

These results have two important implications: i) 
the close correlation observed between S S fr(i.4 ghz) an d 
^sfr(h„) for face-on galaxies reinforces the use of the ra- 
dio luminosity as SFR indicator not only on global but 
also on local scales; ii) extinction could significantly affect 
estimates based on H Q emission for high SFRs in edge-on 
galaxies. 
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Although the star formation rates deduced from the 
H Q emission are systematically underestimated compared 
to those deduced from the radio continuum, the mean SFE 
reported by RY99 for their entire sample (121 galaxies) is 
4.3%, i.e. higher than the mean SFE deduced from the 
radio continuum for our entire sample of 180 galaxies. 
However, considering the 103 galaxies we have in common 
with RY99, the mean SFE deduced from the H Q emission 
is about 3.1% in good agreement with our value. 

5.2. Non-thermal radio continuum scale-length 

Discrete supernova remnants (SNR) themselves account 
only for < 10% of the radio luminosity of normal galax- 
ies (Ilovaisky & Lequeux 1972). Most of cosmic rays es- 
cape from their parent SNRs and diffuse along the galaxy 
disk during their characteristic lifetimes which are much 
larger than typical ages of SNRs (~ 10 5 yr). Therefore 
one expects the non-thermal radio continuum emission to 
be smoothed over a larger region around the parent stel- 
lar population with respect to the H a emission. On scale 
smaller than the cosmic rays characteristic diffusion scale, 
Z?cr, the non-thermal radio continuum cannot provide a 
reliable estimate of the local SFR. On the other hand, 
averaging over regions larger than the typical spiral arm 
width leads to a dilution of the H Q emission. Fig. 12 (top 
panel) shows the ratio between the SFR surface densities 
deduced from H Q and from radio luminosity versus lin- 
ear size of the observing beam. One notes that, apart for 
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Fig. 11. Ratio of SFRs surface densities deduced 
from H„ and 1.4 GHz luminosity for highly in- 
clined (NGC1055 and NGC3079) and starburst galaxies 
(NGC2146, NGC3034). These examples dramatically il- 
lustrate that the SFRs surface density inferred from the 
H Q emission towards the nuclei of these objects are un- 
derestimated by more than one order of magnitude. 
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Fig. 12. Ratio of SFRs surface densities deduced from 1.4 
GHz luminosity and H Q emission as functions of the linear 
scale of the beam (top panel) and distance from galaxy 
center (bottom panel). Highly inclined galaxies (i > 70°) 
are represented with open dots. 



edge-on galaxies (where the H Q is affected by extinction), 
the ratio £sfr(i.4 gh z ) to Xsfr(h„) tends to increase with 
the beam size. In particular, the ratio is smaller or greater 
than one for beam sizes respectively below and above <~ 5 
kpc. This is consistent with a scenario in which the charac- 
teristic linear scale of both cosmic rays diffusion scale and 
arm size is about 3 kpc. Higher resolution radio images 
are required to determine the lowest linear scale at which 
the radio continuum correlates with the H Q emission and 
thus providing an estimate of Dqk- 

We found also that the ratio £sfr(i.4 GHz) to Sg FR ( Ho ) 
on average does not depend on the distance from the 
galaxy centers (Fig. 12, bottom panel). 

5.3. The nature of Schmidt law scatter 

The dispersion (standard deviation) of the composite ra- 
dio Schmidt law shown in Fig. 9 is a factor of 2.8, which 
is smaller than the corresponding dispersion of a factor 4 
of the H Q Schmidt law (RY99). This is probably because 
the radio continuum does not suffer the effects of dust 
extinction. The noisiness of the correlation could reflect 
real variations of the mean Schmidt law, indicating that 
it should be regarded at best as a first order approxima- 
tion of the star formation process, or it could be due to 
the many assumptions involved in the calculation of gas 
and SFR densities. 
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Variations of the Xco factor from galaxy to galaxy can 
be advocated to explain a part of the correlation scatter. 
These could introduce uncertainties up to a factor of 2 in 
the gas density scale (Kennicutt 1998). Another possibility 
is that the extrapolation of the proportionality between 
CO luminosity and virial mass of giant molecular clouds 
observed in our own Galaxy and in nearby galaxies, which 
is the basis of molecular mass determinations in this and 
similar works, does not hold for all spiral galaxies. 

However, by comparing SFR surface densities deduced 
from H Q and from radio continuum luminosity we showed 
that, even in the face-on subsample, the E SFR ( Hq ) — 
^sfr.(i.4 GHz) relation itself is affected by a scatter of at 
least a factor of 2. Hence, the uncertainties of the SFR 
indicators alone can account for a consistent fraction of 
the Schmidt law scatter. 

Fig. 13 shows the Schmidt law separately for Sb, Sbc 
and Sc and SO-Sab, Scd and Sd-Sm galaxies (only detec- 
tions with signal-to-noise ratio greater than 2a are con- 
sidered). The two subsamples have the same mean SFE of 
3.5% but very different dispersions of a factor of 2.4 for the 
former and 3.8 for the latter, i.e. the scatter of the Schmidt 
law is considerably reduced excluding extreme morpholog- 
ical types. RY99 also found that the H Q Schmidt law is 
considerably tightened with the exclusion of the irregu- 
lar galaxies and mergers. These SFE variations around 
the mean Schmidt law can be attributed to the particular 
physical conditions and/or environments experienced by 
these objects (see e.g. starburst galaxies), but they could 
also be induced by observational effects further amplified 
by poor statistic. In fact, Scd galaxies which are char- 
acterized by the lower mean SFE in our sample (see also 
Fig. 5) behave consistently with other morphological types 
excluding IC342 for which the NVSS misses flux density 
from the extended structure, see Sec. 2. 

5.4. Molecular gas consumption timescale 

The star formation efficiency can be expressed in terms 
of gas consumption (or cycling) timescale. The molecular 
gas consumption timescale, th 2 , indicates the time needed 
to convert all the molecular gas in stars given a constant 
SFE. Expressed in Gyr, th 2 is given by 



m 2 (Gyr) 



10 



SFE r8 (%) 



(8) 



The mean star formation efficiency found in this work 
for spiral galaxies yields a molecular gas consumption 
timescale of about th 2 — 2.8 Gyr. Given the dispersion 
of the SFE in the sample, most galaxies have gas cycling 
timescale between about 1 and 9 Gyr. At the lower end 
of the distribution we found galaxies with a SFE Tg ^ 1%, 
such as NGC4535 (see Fig. 2). The corresponding molecu- 
lar gas consumption timescale is almost comparable with 
the Hubble time. The shortest gas cycling timescale cor- 
responds to galaxies with exceptional SFE like NGC3310 
(see Fig. 3). Star formation rates in these galaxies are so 
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Fig. 13. Radio Schmidt law for Sa, Sbc and Sc (top panel) 
and SO-Sab, Scd, Sd-Sm and merging-irregular (bottom 
panel); limits are not shown. Dashed and dotted refer- 
ence lines indicate the mean and the dispersion of the two 
sub-samples, respectively. The scatter of the Schmidt law 
is considerably reduced excluding extreme morphological 
types. The scatter of the correlations shown in top and 
bottom panels is a factor of 2.4 and 3.8, respectively. In 
top panel, the deviation from the correlation of the four 
points belonging to NGC1068 is due to the AGN-related 
radio emission of this Seyfert galaxy (see Sect. 4.2). 



high, compared to their gas supplies, that the gas cycling 
timescale is th 2 < 0.1 Gyr. 

It is interesting to investigate the behaviour of the so- 
called starburst galaxies with respect to the gas cycling 
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Fig. 14. Esfr — Sh 2 plane for all the detection in the 
sample (limits are not shown) along with three reference 
lines corresponding to different molecular gas consump- 
tion timescales, as indicated. Pointings of known starburst 
galaxies are represented by open dots. 



timescale. In literature, galaxies have been classified as 
starbursts according to different criteria. Heckman et al. 
(1998) define a galaxy as starburst when it is hosting a 
star-forming event that dominates its bolometric lumi- 
nosity, i.e. on the basis of the magnitude of the SFR. 
Alternatively, Shu (1987) and Young (1987) proposed a 
classification based on the efficiency of the star formation. 
In this latter definition a galaxy with a high SFR is not de- 
fined as starburst if the mass of gas available is enough to 
sustain the star formation rate and vice versa. Following 
RY99 we selected the galaxies hosting a region in which 
the SFE is enhanced by a factor of three compared to the 
mean of the sample; for these starburst regions th 2 < 1 
Gyr. Fig. 14 shows the usual Esfr — £h 2 plane for all the 
detection in the sample along with three reference lines 
indicating the gas cycling timescales th 2 =0.1, 1 and 10 
Gyr. Most galaxies known in literature as "starbursts" 
have consumption timescales comparable with those of 
normal spiral galaxies. In some cases, e.g. M82, starburst 
galaxies host both regions characterized by SFE lower and 
higher than the mean of the sample. On the other hand, 
galaxies with normal or even low SFR surface densities, 
such as the highly inclined galaxy NGC4631, should be 
regarded as starburst according to our selection criterion 
based on the SFE. Finally, there are some galaxies for 
which the SFE is so high that the gas cycling timescale 
is t H2 < 0.1 Gyr, e.g. the starburst NGC3310. NGC3310 
is one of the best examples of a local UV-bright starburst 



(see Conselice et al. 2000 and reference therein). The ex- 
ceptional high SFE of this galaxy indicates that its inter- 
stellar medium is truly affected by an extraordinary star 
formation event. 

5.5. Conclusions 

The results of this work are the following: 

1. There is a tight correlation between the 20 cm non- 
thermal radio continuum and the CO line intensity in a 
representative sample of 180 spiral galaxies. The correla- 
tion holds within and among the galaxies. 

2. The mean star formation efficiency, i.e. the ratio 
of the radio SFR to the molecular gas densities, for our 
sample is 0.035 x 10 -8 yr _1 with a dispersion of a factor 
of 3. This corresponds to convert 3.5% of the available gas 
into stars on a time of 10 8 yr. 

3. The comparison of SFRs surface densities deduced 
from 1.4 GHz luminosity and from the H a emission for 
102 galaxies, reveals that £sfr(i.4 GHz) and S SFR (h q ) are 
closely correlated for face-on galaxies (i < 40°), reinforc- 
ing the use of the radio radio luminosity as SFR indicator 
not only on global but also on local scales. SFRs surface 
densities deduced by the H Q luminosity for highly inclined 
galaxies (i > 70°) are systematically underestimated for 
Ssfr(i.4 GHz) > 3 x 10~ 3 M Q yr^kpc- 2 . 

4. The star formation efficiency varies weakly (less than 
25%) with the Hubble morphological type. 

5. The variation of the SFE within individual galaxy 
disks is less than a factor of 3. The largest variations are 
found in starburst galaxies. 

6. The SFE is found to be approximately constant as 
a function of distance from the galaxy centers. 

7. The composite radio Schmidt law, star formation 
versus molecular gas content, extends for more than 3 or- 
der of magnitude with an exponent of 1.3. 

8. Most galaxies known in literature as "starbursts" 
have consumption timescales comparable with those of 
normal spiral galaxies. In some cases, e.g. M82, starburst 
galaxies host both regions characterized by a SFE lower 
and higher than the mean of the sample. Furthermore, 
there are some galaxies for which the SFE is so high that 
the gas cycling timescale is th 2 < 0.1 Gyr, e.g. NGC3310. 
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